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Abstract 

We discuss 0{N) invariant scalar field theories in 0+1 and 1 + 1 space-time 
dimensions. Combining ordinary "Large N" saddle point techniques and simple 
properties of the diagonal resolvent of one dimensional Schrodinger operators we 
find exact non-trivial (space dependent) solutions to the saddle point equations 
of these models in addition to the saddle point describing the ground state of 
the theory. We interpret these novel saddle points as collective 0{N) singlet 
excitations of the field theory, each embracing a host of finer quantum states 
arranged in O(N) multiplets, in an analogous manner to the band structure 
of molecular spectra. We comment on the relation of our results to the classical 
work of Dashen, Hasslacher and Neveu and to a previous analysis of bound states 
in the 0(N) model by Abbott. 
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I Introduction 



Field theories involving a very large number N of components have turned out to 
be an extremely useful tool in addressing many non-perturbative aspects of quantum 
field theory and statistical mechanics [[[], 0, |3|. In such theories the action is pro- 
portional to iV (or in some cases, to iV 2 ) which plays the role of l/h. Therefore, as 
N — > oo, the path integral of such field theories is dominated by the saddle points of 
the action, which are protected against being washed out by quantum fluctuations in 
this limit. 

An important subclass of such field theories are the O(N) vector models whose 
dynamical variables are N real scalar fields arranged into an O(N) vector <fr with 
O(N) invariant self interactions of the form (<J? 2 ) n [] which have been studied inten- 
sively in the past 0, ||, || [K|. More recently, 0(N) vector models have attracted 
interest 0, [13|, [14|, |TJ, [16|] in relation to the so called "double scaling limit" || . The 



cut-off dependent effective field theory has also been discussed very recently |TJJ. The 
main object of these studies were the true vacuum state of the O(N) model, and the 
small quantum fluctuations around it. In the case of quartic interactions, of interest 
to us here, the theory is best analyzed by introducing into the lagrangian an auxiliary 
field a that can be integrated out by its equation of motion a= ^ <fr 2 (where g 2 /N 
is the quartic coupling in the lagrangian) .0 In this way the original fields $ may 
be integrated out exactly, yielding an effective non-local action S e ff[a] in terms of o 
alone. 

The vacuum and the low lying states above it are described by the small fluctu- 
ations of the u field around a certain local minimum of S'e/zblQ- in this minimum 
the a field configuration is homogeneous (i.e. space-time independent) a = a c , where 
the constant o c (which is the true mass of the $ "mesons" around that vacuum) 

1 We concentrate on the case of quartic interactions n = 2 

2 For higher interactions we need an additional lagrange multiplier field. 

3 This minimum may not be the absolute minimum of S e ff[a] in certain circumstances flTo| , |jj 
but the decay rate of the metastable vacuum is suppressed by a very small tunnelling factor of the 
order of e~ NV , where V is the volume of space-time. 



1 



is obtained from the gap equation — \ a=a = 0. Homogeneity of the vacuum a 
configuration is expected due to Poincare invariance of the latter. 

S e ff[a] is a rather complicated functional of a and in addition to o = o c it possesses 
yet many other extremal points at which the solution to the extremum condition 
5Se sl ^ = are space-time dependent a configurations a(x). These extrema cannot 
describe the ground state (and indeed, neither a metastable ground state) and are 
thus saddle points of S e ff[o~] rather than (local) minima. Indeed, it can be shown by 
direct calculation^ that the second variation around them has an infinite number of 
negative eigenvalues. They correspond to collective "heavy" O(N) singlet excitations 
of the <E» field, where its length squared oscillates in isotopic space as a function of 
space and time. 

In this work we find explicit exact formulae for these o~(x) configurations in + 1 
and in 1 + 1 space time dimensions. In the two dimensional case we discuss only time 
independent a configurations. In either dimensionalities they have the general form 
a(x) = ip(x — Xo), where ip is an elliptic function and xq is an arbitrary parameter 
responsible for translational invariance. In + 1 dimensions if the amplitude of ip(t) 
along the real axis is finite it corresponds to a finite S e ff[cr] (after regularizing the 
contribution of the fluctuations of the O(N) vector field into a finite expression), 
while if the amplitude is infinite (i.e. the pole of ip lies on the real axis) it is a 
potential source of an infinite S e ff[a], which can suppress this mode enormously 
relative to the homogeneous a = a c configuration. Moreover, such infinite amplitude 
a(t) configurations, imply very large fields, a fact that might throw us away from 
the validity domain of the large N approximation. We will always keep this in mind 
when addressing such configurations, but in + 1 dimensions we argue that they 
might be relevant in certain circumstances. In the two dimensional field theory case 
most a(x) = ip(x — x ) configurations give rise to infinite S e ff[a] values (whose 
divergences are inherent and cannot be removed by regularization) and are therefore 
highly suppressed. Infinite amplitude o~(x) configurations fall obviously into this 
4 This will not be carried out here. 
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category, but so do finite amplitude o~(x) configurations. Some of the latter yield 
(regularized) finite S e ff[a] values only in the limit in which their real period (as 
elliptic functions) becomes infinite. Similar o~(x) configurations were found in the two 
dimensional O(N) model in |22| by using inverse scattering techniques. Nevertheless, 
they turn out to be quite distinct from those discussed here. Our finite amplitude 
periodic o~(x) configurations might become important if we put the field theory in a 
spatial "box" of a finite length L, but as L — > oo they will be hightly suppressed due 
to their diverging S e ff[a]. 

These a{x) saddle point configurations are hybrid objects, in the sense that they 
are static solutions of the classical equations of motion SSe g* — = of the action 
S e ff[a] that embodies all quantal effects of the $ field. Therefore, due to all their 
properties enumerated above, the finite S e ff[a] a(x) configurations are reminiscent, 
in some sense, of soliton solutions in other field theories, but this analogy is only 
a remote one, since the <j(x)'s carry no topological charge of any sort and they are 
not local minima of S e ff[o~]. They are more similar to the so called "non-topological 



solitons" |24j, p5| . In this paper we ignore the gaussian fluctuations of the a field 
around the saddle points cr(x), but explain how to find the spectrum of "mesons" (i.e. 
$ quanta) coupled to the o~(x) configuration. It is clear that each o~(x) configuration 
gives rise to a tower of "meson" bound states as well as to "meson" scattering states, 
arranged in highly reducible O(N) multiplets. The infinite negative eigenvlues of 
the second variation of S e ff[o~] around these configurations correspond to the infinite 
number of decay modes of such excited states into all lower ones. This form of 
the spectrum of quantum states is reminiscent of the vibrational-rotational band 
structure of molecular spectra |TJ. Since we ignore fluctuations of the a field around 
the saddle point configuration o~(x), the "molecular analog" to this would be to have 
the molecule only at one of its vibrational ground states and consider the rotational 
excitations around it. As a matter of fact, in the + 1 dimensional case, which is 
nothing but the quantum mechanics of a single unharmonic oscillator in N euclidean 
dimensions this analogy is actually the correct physical picture. In this case the o~(t) 
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configurations are actual vibrations of the iV dimensional vector, changing its length. 
As the vector stretches or contracts, it might rotate as well, changing its orientation 
in N dimensional space, which is energetically much cheaper. The latter gives rise to 
the rich O(N) structure of the spectrum. 

Our approach to the specific problem of bound state spectrum in the 1 + 1 dimen- 
sional O(N) model has been inspired by the seminal papers of Dashen, Hasslacher and 
Neveu (hereafter abbreviated as DHN) fl8fl , [|H| , p0| , and especially, by their analy- 
sis of the bound state spectrum in the Gross- Neveu model In [^(]] DHN present 
a detailed construction of a sector of the bound state spectrum of the Gross-Neveu 
model whose corresponding u{x) configurations^ differ only slightly from the homoge- 
neous vacuum configuration o = <j c of the Gross-Neveu model, at least as far as states 
with a principal quantum number well below iV are concerned. For such states, the 
corresponding a(x) configuration looks like a pair of interacting (very close) kink and 
anti-kink, where the kink amplitude is of the order 1/N. This might be thought of as 
the back-reaction of the fermions in the bound state on the vacuum a = a c polarizing 
it into a a[x) configuration. Despite the smallness of a(x) — cr c , it has a dramatic 
consequence: the production of its corresponding bound state. Indeed, DHN find the 
explicit form of cr{x) by performing an inverse scattering analysis of the one dimen- 
sional Schrodinger operator d 2 + a(x) 2 — a 2 + <j'(x) obtained straightforwardly from 
the Dirac operator %$ — a. The potential term in this Schrodinger operator vanishes 
identically and cannot bind for a = a c , but for the a(x) configuration of DHN it 
becomes an attractive reflectionless potential (of depth of the order 1/N 2 ). Thus, 
despite the small magnitude of the "back-reaction" of the fermions on the vacuum, 
it is enough in order to create a very shallow dip in the one dimensional Schrodinger 
potential, causing it to produce one bound state. This is the reason why DHN at- 
tribute the appearance of fermionic bound states in the Gross-Neveu model to its 
intrinsic infra-red instabilities (the distortion of a c into cx(x)) which result from the 
asymptotic freedom of the model. Indeed, the same infra-red instabilities generate 

Here <j(x) stands for the iptp bilinear fermionic condensate. We refer here only to their static 
<j{x) configurations. 
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non-perturbatively the fermion mass ga c in this model in the first place. 

The 1 + 1 O(N) vector model with quartic interactions being trivially "asymp- 
totically free" due to its super-renomalizability, exhibits for negative quartic coupling 
a sector of bound states (whose a(x) configuration was mentioned above) that is very 
similar to the bound states of the Gross-Neveu model. It is probably created by the 
same mechanism of vacuum infra-red instabilities. Their energy spectrum and corre- 



sponding reflectionless kink like a(x) configurations were found in |22 by following 
the DHN prescription step by step. 

In addition to the description above, DHN mention briefly yet another sector of 



bound states [2(| in the Gross-Neveu model, namely, those built around the Callan- 
Coleman-Gross-Zee kink. The latter is a local mimimum cr(x) configuration very 
different from the homogeneous vacuum one. It connects the two true degenerate 
vacua of the Gross-Neveu model. Here all fermions are bound at zero binding energy 
in the fermionic zero mode of the kink (which is the only bound state supported by 
the kink), with no back-reaction at all on the kink independently of the number of 
fermions trapped in that state. The cr(x) configurations we find in the O(N) model 
are precisely of this type. However, the latter cannot be associated with degenerate 
vacua, since the O(N) model lack such a structure. Due to this fact and since we 
do not use inverse scattering techniques, our results are complementary to those of 



20, 22]. Our method could be used also to find analogues to the Callan-Coleman- 



Gross-Zee kink in other two dimensional theories, where inverse scattering methods 
become practically useless [p8 |. As far as we know this is the first use made of 
the method described below. The idea we use is very simple. A generic saddle point 

SS [crl 

condition for the O(N) model (in any dimension) — e ^ = relates a to the diagonal 
resolvent of a Klein-Gordon operator (for scalar fields) with a as its potential term. 
In the static case the Klein-Gordon operator reduces into a Schrodinger operator. In 
1 + 1 dimensions the resulting diagonal resolvent of the one dimensional Schrodinger 
operator is known to obey a certain differential equation, known as the Gelfand-Dikii 
equation [2~3|. Since the resolvent is essentially cr(x) due to the saddle point condition, 
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the latter must obey an ordinary differential equation induced from the Gelfand-Dikii 
equation. 

The differential equation imposed on a(x) is of second order. Solving it for cr(x) we 
encounter two integration constants. One constant is trivial and insures translational 
invariance of the solution for <j(x) while the other acts as a parameter that distorts the 
shape of a(x) and parametrises its period along the x axis. When cr(x) has a finite 
amplitude we can always use this parameter to drive u{x) into an infinite period 
behaviour. In this limit <j(x) reproduces the reflectionless configuration analogous to 
the Callan-Coleman-Gross-Zee kink in the case of negative quartic coupling. 

At this point it is worth mentioning that some bound states in O(N) models in 
various dimensions were found around the homogeneous vacuum a = a c configuration 
in the singlet $ • $ channel of "meson-meson" scattering, by simply analyzing the 
poles of the a — a propagator 0, 0, [|Kj.[] Specific conditions for the existence 
of such bound states and resonances in two space-time dimensions to four may be 
found in || [|. In certain theories in three space-time dimensions such bound states 
may become massless (to leading order in 1/N) and are identified as goldstone poles 
(a dilaton 5fJ or a dilaton and a dilatino P?[ ) associated with spontaneous scale 
invariance breakdown. Occurance of such massless bound states is responsible for the 
double scaling limit in vector models fl6| and the consequences of obstructions for 
this to happen in two space-time dimensions were addressed in ||14 . 



The paper is organized as follows: In section (II) we analyze the + 1 dimensional 
model. This is done mainly to introduce notations and gain confidence in our method. 
In particular, we show that our method for obtaining the time dependent saddle point 
a configuration is completely equivalent to a JWKB analysis of the radial Schrodinger 
equation of the unharmonic oscillator where h — 1/N establishing the validity of our 
method. We analyze both cases of real and imaginary time qualitatively, without 
giving explicit formulae for the a(t) configurations. 

In section (III) we turn to the 1 + 1 dimensional case. Here as well we analyze both 
6 See also the third reference in 0. 
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Minkowsky and Euclidean signatures of space-time. Our static a(x) configurations 
are obtained in a certain adiabatic approximation to the saddle point equation. Non- 
zero frequency modes are strongly suppressed and the only meaningful result of the 
saddle point equation concerns the zero frequency mode. In this way we effectively 
reduce the theory from 1 + 1 dimensions to + 1 dimensions. We give an explicit 
expression for some of the cr(x) configuration that is important in the infinite volume 
case and construct the spectrum of "mesonic" bound states in its background. 

We draw our conclusions in section (IV). A simple proof of the Gelfand-Dikii 
equation, emphasizing its elementarity in the theory of one dimensional Schrodinger 
operators is given in an appendix. 
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II The Unharmonic Oscillator 



We begin our investigation of the bound state spectrum in the two dimensional 
O(N) vector model by considering the simpler case of the quantum mechanics of 
a single O(N) invariant unharmonic oscillator. Clearly, the most straightforward 
way of analyzing the quantum state spectrum of this system when iV — ► oo is to 
apply the large N approximation directly to the hamiltonian of this oscillator. This 
procedure is entirely equivalent to the JWKB approximation to the corresponding 
radial Schrodinger equation, where % = ■h Q . Comparing our method of calculating 
the spectrum, which is an analysis of the action in the large N limit, to the JWKB 
approximated Schrodinger equation reveals their equivalence, which is not surprising, 
since the large N approximation to the path integral is equivalent to the semiclassical 
approximation with the role of h played by 1/N. 

Therefore, our purpose in this section is to better understand the way our method 
works in this simpler case and gain confidence in it before generalizing to the field 
theoretic case. Such an understanding is achieved precisely by demonstrating its 
equivalence to the JWKB approximated Schrodinger equation. This reasoning is 
completely analogous to the one made in [^, [j] , in order to understand the role of 
instantons as tunnelling configurations in the path integral. Indeed, tunnelling effects 
in quantum mechanics are calculated straightforwardly by using JWKB approximated 
Schrodinger equations rather than path integral instanton methods, but it is the 
instanton calculus that is generalizable to quantum field theory and not the JWKB 
approximated Schrodinger equation. Since our main interest in this paper is the 1 + 1 
dimensional field theoretic case (discussed in the next section) we will demonstrate the 
consistency of our method in the + 1 dimensional case without giving the resulting 
a(t) configurations explicitly in terms of elliptic functions of t. Evaluation of these 
functions may be done straightforwardly from our discussion. 
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The O(N) invariant unharmonic oscillator is described by the action[] 

T 

S = [ dt \lx 2 - a (^m 2 x 2 + ^ (3g 2 (x 2 ) 2 



(2.1) 



where x is an N dimensional vector, g 2 is the quartic coupling and m is the pure 
harmonic frequency (as g = 0) of the oscillator. Here a and j3 are discrete parameters 
whose values are either +1 or — 1. a — 1 corresponds to the real time case, while 
a = — 1 corresponds to imaginary (euclidean) time. Likewise, (5=1 corresponds 
to the well defined stable case where the quartic potential goes to +oo as \x\ — > 
oo, while (3 = — 1 corresponds to the upside down potential. We have introduced 
these parameters for efficiency reasons, allowing us to discuss all different possibilities 
compactly. In the large iV approximation g 2 is finite, hence the quartic interactions 
in Eq. ( |2.1| ) are in the weak coupling regime^. 
Rescaling 



X' 



Ny< 



(2.2) 



Eq. fl2.ip may be rewritten as 



S 



T 

N J dt 
o 



1 



y 



m 



a 



2\2 



(y 2 ) 



(2.3) 



Introducing an auxiliary variable a, the action is finally recast into 

'1 . 



S 



T 

N f dt 



o 



am 



2 V 



y +a(3(a -gay' 



(2.4) 



which is the starting point of our calculations. Note that a has no kinetic energy and 
thus may be eliminated from Eq. ( f2.4[ ) via its equation of motion 

o = \gy 2 - (2.5) 



We assume throughout this work that m 2 > 0. A sign flip of m 2 is equivalent to a simultaneous 
sign flip of the parameters a and f3. 

8 Taking the scale dimension of m to be canonically f the other dimensions are [t] = — 1 , [x] = 

~\i L? 2 ] = 3- The dimensionless ratio 9 Jj? which is much smaller than 1 for finite g 2 ,m 2 clearly 
implies weak coupling. 
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Throughout this paper we adopt the convention g > and therefore a is a non- 
negative variable. 

Since in principle we are interested in calculating the spectrum of bound state 
of the oscillator we consider the partition function Tre~ HT = J d N y Q Wp(T) where 
Wp(T) is the transition amplitude of the oscillator to start at y and return there 
after time-lapse T 

(y ;t = T\y ; t = 0) = W P (T) = J VyVa e^ s . (2.6) 

pbc 

Here the subscript u pbc" denotes integration over paths obeying the obvious periodic 
boundary conditions implied by Eq. ( |2.6| ) and 7 = % for the real time case, while 
7 = — 1 for imaginary time. Clearly, the periodic boundary conditions on o are 
dictated by Eq. ([T5D. 

Eq. ( p.4|) is quadratic in y which may be therefore integrated completely out of 
the action leading to 

W P (T) = J " Va det;^ /2 [-d 2 t - am 2 - 2a(3ga] e N ^Io^ dt (2.7) 

pbc 

where the operator whose determinant is taken in Eq. ( j2.7|) is defined on the interval 
[0, T] with periodic boundary conditions. Exponentiating the determinant in Eq. 
( |2.7| ) we may express Wp(T) as 

W P (T) = J Va e N ^ff [(T] (2.8) 

pbc 

where 

S e ff[cr] = a(3 / a 2 — — Tr p b c ln — d 2 — am 2 — 2a(3ga (2.9) 
J 27 L - 1 

' 

is the exact effective non-local action in terms of the a field. 

Due to the explicit factor of iV in the exponent of Eq. (|2.8j), Wp(T) is dominated 
by the extremal points of S e ff[a], and to leading order in 1/N f\ it is given by the 
9 Thus ignoring quadratic fluctuations of a around these extrema. 
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coherent sum over all such extremal configurations o~(t) of the integrand in Eq. ( |2.8| ) 
evaluated at these configurations pQ| . 

These extremum configurations are solutions of the functional condition 

6S eff [a]/S(7 = ; a(0) = a(T) (2.10) 

which is equivalent to 



g_ 

7 x 1 —d 2 — am 2 — 2a(3ga 



2a(t) = - y -(t\ \t); (7(0) = a{T) . (2.11) 



Eq. (|2.11|) is a rather complicated functional equation for generic a configurations. 
However, it simplifies enormously for constant a = o c configurations which obey the 
periodicity boundary condition automatically for any T. In this case the diagonal 
matrix element is given by a simple Fourier integral (for T — > oo), and Eq. (|2.11|) 
descends into the cubic equation 

S2(3go* + 16m 2 a 2 - g 2 = . (2.12) 

Eq. ( |2.12j ) has been used in JTT, [TJ|, PI to determine the optimal frequency squared 
uj 2 = m 2 + 2j3ga c of the harmonic approximation to Eq. (|2.1|) for low lying excitations 
around the vacuum (or metastable vacuum, for (3 = —1) of the unharmonic oscillator. 

Following Eq. ( |2.5| ) one should obviously consider only positive real roots of Eq. 
(|2.12|) . Analyzing the left hand side of Eq. ( |2.12|) it is clear that for /3 — +1 there 



is only one such root, while for (3 = — 1 there are two, but only the smaller one 
corresponds to the metastable vacuum.0 These assertions will be demonstrated more 
clearly when we will examine the JWKB approximation to the radial Schrodinger 
equation of the unharmonic oscillator in the last part of this section. 



10 



In case of finite dimensional integrals, not all saddle points necessarily contribute. Only those 



that are located along the steepest contour contribute [Bl[. In the functional integral (Eq. (|2.8|)) 
the situation is not so clear. However all t dependent extrema a(t) we find have simple classical 
meaning, and therefore must be all important. 
11 See also the third reference in 0. 

12 In this case uj 2 c = m? — 2ga c and uj 2 < for the larger root. 
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We now turn to generic time dependent solutions of Eq. ( |2.11|) . This equation 



relates the diagonal resolvent R(t) = (t\(h — am 2 ) 1 \t) of the Schrodinger operator 

h = -d 2 - 2a(3ga (2.13) 

whose potential term is U(t) = —2aj3ga, to o~(t). This diagonal resolvent obeys the 
Gelfand-Dikii equation (Eq. (A. 12)) discussed in the appendix. Thus, eqs. (|2.11|) 
and (A. 12) imply 

2aa" - (a 1 ) 2 + Aa(m 2 + 2f3gaW = ^- ; a(0) = a(T) (2.14) 

4a 

where a' = da/dt and we have used the relation a = — 7 2 . 

Eq. ( |2.14| ) is one of the main results of this paper. It implies that the functional 
equation Eq. ( 2.11| ) is equivalent to an ordinary differential equation that we readily 



solve below, obtaining the required extremal a(t) configuration. 

Note that Eq. (|2.14j) admits a time independent solution a = a c for which it 
reproduces Eq. ( |2.12| ) exactly. Regarding time dependent solutions, first integration 
of Eq. (|2.14|) may be done straightforwardly by substituting 

a'(t) = f(a) (2.15) 



which leads to 



9 2 (o_ 
0o 



4 (m 2 a 2 + (3ga^ 



<r(0) = a{T) (2.16) 



where ctq is an integration constant (and we have used a 2 = 1). Eq. ( |2.16| ) has the 



form of the equation of motion of a point particle in one dimension along the ray[^] 
cr > in the potential V(a) given by 

- V(a) = 9 - - l) - 4 (mV + Pga 3 ) ; a > (2.17) 



13 H ere w e assume that the constraint a > stemming from the classical equaiton of motion of a 
(Eq. (|2.5|) remains in tact quantum mechanically. We will see that this is indeed the case. Since 
Eq. ( p.5[) could be enforced as an exact functional constraint by introducing a lagrange multiplier 
field, this conclusion should not be of any surprise. 
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(we have shifted the constant piece of V such that Eq. (|2.16| ) corresponds to a motion 
at energy zero). Note that a has been factored out in front of — V in Eq. fl2.16|) , as it 
should, since a = +1 corresponds to motion in real time while a = — 1 corresponds 
to motion in imaginary time. The classical turning points are the non-negative roots 
of the cubic equation 

V(a) = 0. (2.18) 

For real time motions, the classical a paths will be confined to regions where 
— V > 0, while for imaginary time motions, to regions where —V < 0. In addition to 
the turning points given by Eq. ( p.l8j ), there is a turning point at a = due to the 
constraint o > 0.0 

Denoting the roots of Eq. fl2.18p by 01, cr 2 , cr 3 Q the following formulae are useful 
for listing the various solutions of Eq. ( |2.16|) . 

a 2 

-1/(0) = -j<0 (2.19a) 
-V(±oo) = TP-oo (2.19b) 
o"io"20"3 = ~y| • (2.19c) 

Assuming m and g are given, o"o is the only arbitrary parameter, arising as an inte- 
gration constant of Eq. ( |2.16| ). We now show that equation (|2.16|) is equivalent to the 
(one dimensional) radial equation of motion of the oscillator in a quantum effective 
potential obtained as N — > oo, whose energy parameter £ is inversely proportional 

tO 0Q. 

This observation is conceptually an important one, assuring the validity of our 
calculation of the time dependent extremal a configurations of S e ff[a], identifying 
them as corresponding to collective radial excitations of the oscillator, as we have 
discussed at the beginning of this section. Moreover, holding m and g fixed, cri,cr 2 
and o"3 are functions of o"o, and in the most generic case two of them are complex 



14 cr trajectories that reach the endpoint a — always occur in imaginary time. They are unde- 
sirable since they have no quantal counterpart. We comment on this point in the last part of this 
section. 

15 We adopt the convention that when all three roots are real they are ordered as u\ < 02 < 0-3. 
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(conjugate). Therefore, the best way to analyze the various turning point configu- 
rations of the solutions to Eq. (|2.16|) is to consider the (large N) effective radial 
potential of the oscillator giving rise to these radial motions. Varying the energy 
parameter of the latter will show us the way the turning points of Eq. ( |2.16 ) change 
as functions of ctq. Recall that Wp(T) defined in Eq. ( |2.6| ) is given as a sum over all 
saddle points a(t) found from Eq. ( [2.11| ). These are essentially parametrized by a$. 
Thus a crude estimate of Wp(T) would be to insert these saddle point configurations 
into the integrand of Eq. ( |2.8| ) and sum over all allowed values of do- Since <To turns 
out to be inversely proportional to the energy parameter £ of the associated radial 
Schrodinger equation this is just the natural thing to do, since by definition Wp(T) is 
a Laplace transform of the transition amplitude Wp(S) in the energy plane. In what 
follows, howerver, we will not calculate Wp(T) but rather concentrate on a specific 
saddle point a(t) with its particular cr parameter. 

Before turning to the hamiltonian formulation we depict qualitatively in figures (1) 
and (2) the turning point structure of Eq. ( [2.16|) in the case where o~i, o 2 and 03 are all 
real. Figure (1) corresponds to the case (3 = +1 while Figure (2) describes the various 
possibilities for (3 = -1. When f3 = +1, Eqs. imply -V(0) < 0, - V(=F 00) = 

±00 and o\Oio?, < 0. Therefore, if <Ti, <j 2 and 03 are all real, either o\ < a 2 < 03 < 
or o~i < < o"2 < o"3- These two possibilities are shown in figures (la) and (lb), 
respectively. For f3 = — 1 we have — V(0) < 0, — V(±oc) = ±00 and aio^o^ > 0. 
Thus, if 0i, cr 2 and a 3 are all real either <j\ < a 2 < < 03 or < <J\ < a 2 < (t 3 . These 
two possibilities appear in figures (2a) and (2b), respectively. 
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FIGURE 1: The (3 = 1 case of Eq. flTiq) . Either ^ < tr 2 < a 3 < (a) or 
G\ < < ct 2 < <7 3 (b). Bold faced dots correspond to turning points of the a(t) 
motions. The empty circle denotes the reflecting infinite potential wall at o = 0. The 
regions along the a axis accessible to real time motions are marked by the zig-zag 
line and those accessible to motions in imaginary time-by double lines. 




FIGURE 2: The (3 = -1 case of Eq. fl2~lp . Either o x < cr 2 < < tr 3 (a) or 
< <7i < o"2 < 0-3 (b). For various notations see caption to Fig. (1). 
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The Schrodinger hamiltonian corresponding to the action of Eq. ( |2.1|) is 

h 2 



H 



A+U(r) 



(2.20) 



where r 2 = x 2 , A is the N dimensional laplacian and 



U(r) 



m 



— r H 

2 4iV 



(2.21) 



Separating out the radial Schrodinger equation for a state carrying angular mo- 
mentum i we find 



2 r ^2 



«9 2 N-l d i(£ + N-2) 



dr 2 



r dr 



+ U{r)\ R n/ (r) = £ n/ R n/ {r) 



i — >0 



(2.22) 



where n is the principal quantum number. Following j|] we perform the similarity 
transformation 



H^H = r- v Hr v , R n/ {r) - w (r) = r^i^r), i/ = 
which disposes of the term in Eq. (12.221) .^ We finally obtain 

H = Nh = N 

where p 2 = y 2 = r 2 /N and 

V(p) = 
wnAp) = h 2 

Eqs. (ggg ) and (ggg ) imply 



iV-1 



(2.23) 



•^T^+V(P)+«V(P) 



^ 7// x ^ ™ 2 2 /V 4 



(2.24) 



2N 



i H)H) 

2iV 2 



2 ' 



(2.25) 



(2.26) 



16 This transformation if familiar from elementary quantum mechanics. What it does, is essentially 

N — 1 

swallowing a square root of the radial jacobian (the p~^~ factor) into a wa ve fu nction in an inner 
product: R n i(r) — > Xn,e{r), inducing the corresponding transformation Eq. (2.23) on H. 
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assuring the self adjointness of H. 17 For a later reference let us recall that the singular 
solution Xni to the equation Hxn,i — £n,t Xn,t that is definitely not in the Hilbert space 
at hand, blows up at p = asQ 



XuAp) (2.27) 

The explicit ft dependence of V and wn,i in Eq. ( |2.25| ) is obvious- its origins are 
the centrifugal barrier in Eq. ( [2.22| ) and the similarity transformation we performed 
on the operator H. To leading order in 1/N the only ft dependent term that survives 
there is the ^2 piece in V. This is consistent with Langer's modification of the 
centrifugal term ±£(£ + N-2) -> \(i+^) 2 in the limit N -> oo @. It is therefore 
the leading 1/N quantum contribution to the effective potential. Thus, it must be 
equivalent to the leading 1/N correction to the action from the determinant in Eq. 
( |2.9| ). We show now that this is indeed the case, namely, as we have stated earlier 
the solutions of Eq. ( |2.16| ) correspond exactly to classical motions in the potential 
V(p). From now on we set ft = 1 for convenience in eqs. ( |2.24|) -( p.25|) . Clearly, 
from the point of view of the large N approximation, 1/N will play the role of a 



small redefined Planck constant. Therefore 1/N analysis of Eq.( 2.24 ) is equivalent 
to its JWKB analysis. But this is clearly not the ordinary case considered in a 
JWKB approximation, namely, the potential in Eq.( [2.24|) depends explicitly on the 
small parameter 1/N. While this may complicate the 1/N expansion of the actual 
eigenvalues of ft in Eq.( |2.24|) , as far as the 1/N expansion of wavefunctions away from 
the turning points is concerned, the situation is much simpler. Indeed, substituting 
the expansion 



Xe(p) =exp{iN 



Ln=0 



(2.28) 



into the eigenvalue equation for ft in Eq. (|2.24 ) (assuming the eigenvalue is independent 
of N up to an overall scaling) yields to first sub leading order in 1 /N 

x,(p) = [2(£-v( P ))r 1/4 



"Alternatively, Eq. (2.26) may be derived directly from Eqs. fl2.24| )-(|2.25|). 

18 This can be seen immediately from the fact that the wronskian W(x, x) is necessarily a non 
vanishing constant. 



17 



exp | ±i iV f [2 (£ - V W) ] V (i - -j—j + 0(J V-)j «, j . (2.29) 

Therefore, the classical momentum associated with Eq.( |2.29|) is 



p(p) = [2{£-V{p))} 



1/2 



V 2Np 2 2(8 -V) 
[2 (^-V( P )-^(p))] 1/2 + 0(iV- 2 ). 



+ o(i\r 



(2.30) 



Thus, to leading order in 1/N we may safely throw away all non leading terms in the 
potential of Eq.( 2.24|) (which include all £ dependent terms) obtaining 

y4 a (l) 2 = £ - v (") ■ "< r >="<°> 

( |2.25| ) we find that Eq. ( |2.3ip may be written as 



(2.31) 



Using Eqs. 



a 



da 



4 



9 



a - 1 - 4(mV + Pga 6 ), a{T) = a(Q) 



which is nothing but Eq. ( |2.16| ) upon the identification 

16£ _ 1 
9 ' 



(2.32) 



(2.33) 



Thus, we have proven that our calculation of time dependent saddle point con- 
figurations a{t) (Eqs. (|2.14j) - (|2.16| )) of S e ff[a] is equivalent completely to the time 
dependent classical solutions in the effective potential V (Eqs. ( |2.25| )-( |2~3"T| )). Note 
from Eqs. ( [2.24| )-( p.25| ) that angular momentum effects will show up in S e ff[a] only 
in higher terms of the 1/N expansion. Post factum the equivalence of Eqs. ( |2.1ti| ) 
and ( p.31| ) seems self evident, but we have presented an explicit proof of it. More- 
over, in passing from Eq. (|2.31| ) to (|2.32|) we have used the classical a equation of 
motion (Eq.(|2~5|)). This equation must hold also quantum mechnaically because we 
can enforce it as a constraint by introducing a lagrange multiplier field in an equiva- 
lent formulation. Indeed, note that Eq. ( |2.5| ) used in passing from Eq.( 2.31 ) to ( |2.32| ) 
relates the slow radial semiclassical motions p(t) of y in the effective potential V (that 
already includes quantal corrections from y fluctuations) to the a fields. 
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The equivalence of Eqs. fl2.14| ) and ( |2.29| ) holds for their time independent solu- 
tions as well. Indeed, from Eq. (|2.24|) the extremum condition for V(p) reads 

dV 1 
~dp ~ ~4p~ 2 



+ m p + 



fp 3 







(2.34) 



which is equivalent to Eq. ( |2.12|) , the equation determining constant a = o c extremal 
configurations of S e ff[a}. The ground state of the unharmonic oscillator correspond 
to the local minimum of V(p). For (3 = — 1 this is only a metastable ground state. 
It will decay at a rate suppressed by a factor of e~ N to — oo. V(p) has been drawn 
schematically in figures 3(a) and 3(b) (for f3 = 1 and f3 = — 1, respectively) where we 
have also denoted the various regions along the p axis accessible to motions either in 
real or imaginary times for some energy parameter S. Sq = V(po) in Fig. (3a) (Fig. 
(3b)) corresponds to the absolute (metastable) ground state of the oscillator, thus, 
using Eqs. (|2~5l) , (|2.34|) we have 

1 



9pl ■ 



(2.35) 



For /3 — — 1 Eqs. fl2.12|) , ( [2.34|) have yet another solution, corresponding to the local 
maximum in fig. (3b). 



^(P) 





(b) 



FIGURE 3: Radial motions from Eqs. (ggg ) and (|]3T]) as N -»• oo for (3 = +1 
(a) and (3 = —1 (b). For the various notions see the caption to Fig. (1). 



19 



The consistency of Figs. (la,b) and (3a) and of Figs. (2a,b) and (3b) is apparent. 
Clearly Fig. (lb) corresponds to the case £ > £ in Fig. (3a) while Fig. (la) 
corresponds to £ < £o-Q Similarly, Fig. (2a) corresponds to the cases £ > £ c or 
£ < £q in Fig. (3b) while Fig. (2b) corresponds to £ < £ < £ c . 19 

Having figure (3) at hand it is clear now how the turning point structure of Eq. 
( |2.16|) and figures (1) and (2) change under variations of the parameter do holding m 
and g fixed. Consider a classical trajectory in Fig. (3a) with a very large and positive 
energy parameter £. By virtue of Eq. ( |2.33| ) this corresponds to a very small positive 
Co parameter in Fig. (lb). Clearly in this situation we have 0+ ~ oi <C 03. Reducing 
£ in Fig. (3a) (equivalently, increasing cr ), there are no qualitative changes in Fig. 
(lb). It is clear that as £ decreases, o<i increases while 03 decreases, diminishing the 
amplitude of the real-time motion, and in addition, reducing the value of —V(a) at 
its maximum A. These changes persist until £ hits £0 = V(po)> where 02 and 03 
coincide, shrinking the real time trajectory to a point, corresponding to the fact that 
the oscillator is frozen at its ground state. In this case the common value of 02 and 
03 is the constant configuration a = a c mentioned above. Decreasing £ further, we 
flip from Fig (lb) to Fig. (la), where only imaginary time motions are allowed. A 
similar behaviour occurs in figures (3b) and (2a,b). 

Figures (l)-(3) exhibit imaginary time solutions to Eqs. (|2.16|) and fl2.31|) that hit 
the boundary point = for any value of £ or cr - Such solutions describe situations 
in which the particle may be close to the end-point = during a finite time period. 
This is in obvious contradiction with the small p behavior of the square integrable 



wave functions Xn,e(p) in Eq. ( 2.26 ) which suppress the probability of having the 
oscillator oscillating at ~ p <C 1 completely (for any £) as iV — > 00 .Q'Q We will not 
discuss these solutions further. 

As was mentioned at the beginning of this section, a flip in the sign of m? is 



19 With the exception that (71,02 there may become complex (conjugate). 

20 The same problem arises for a free particle in N dimensions and imaginary time. 

21 It might be though that these classical trajectories are associated with the functions \n,i in Eq. 



( |2.27| ) that are solutions of the Schrodinger equation, but do not belong to the Hilbert space. 
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equivalent to a simultaneous flip of a and (3, i.e. a change from real (imaginary) time 
motions in Fig. (3a) to imaginary (real) time motions in Fig. (3b). Therefore, if 
one starts with m 2 > 0, a flip in the sign of m 2 will not induce a local minimum in 
addition to the one that already exists in Figs. (3a) and (3b) at p = po- Alternatively, 
one can check explicitly, that the local minimum of U{p) in Eq. ( [2.21| ) at p = —m 2 /g 2 



(for m 2 < 0) is distorted away by the l/8p 2 term in in Eq. ( |2.25| ). 

We end this section with some general observations and remarks on the possi- 
ble solutions a{t) to Eq. ( |2.16|) . Following our discussion of the various possible 
trajectories, a straightforward integration of Eq. (|2.16| ) yields the latter as 

a(t) =ip(t -t ) (2.36) 



where ip is an elliptic function |34], |33| and to is an integration constant assuring time 
translational invariance. As should be clear from Figs. (l)-(3) the amplitudes of the 
oscillator are either bounded (the poles of ip in Eq. ( |2.36| ) are away from the real t axis) 
or unbounded (i/j has a double pole on the real t axis). In the latter case ip(t — to) is 
essentially a Weierstrass V function, while in the former case, it may be written most 
simply as a rational expression in terms of the Jacobi sn function. These finite am- 
plitude solutions descend into the expected harmonic oscillations in real (imaginary) 
time as £ — > £ + (£ — * £c~) in Figs. (3a) and (3b), where the oscillator executes 
small oscillations around the relevant minimum. The unbounded oscillations on the 
other hand are potential sources to an infinite S e ff[a], which can suppress such modes 
completely. Clearly, as they stand, they are not even periodic in time, since there 
is nothing to reflect a(t) back from infinity. Moreover, these modes are associated 
with very large values of x, which might throw us out of the validity domain of the 
large N approximation. However, regardless of all these problems, the Weierstrass 
V function a(t) configurations are very interesting since for such cr's the Schrodinger 
hamiltonian h in Eq. fl2.13Q , which is essentially the inverse propagator of the y field 
in Eq. ( |2.4Q , is one of the completely integrable one dimensional hamiltonians of the 
Calogero type p5 . 



For real time motions, such trajectories occur only for f3 — — 1 (see Fig. (3)), 
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i.e. only in the case of negative quartic coupling. In such a case the system, strictly 
speaking, does not have a ground state and is meaningless beyond large N saddle 
point considerations. However, due to the very fast decrease of the potential to — oo, 
the time of flight to infinity is finite, and the unharmonic oscillator attains in this case 



a one parameter self adjoint extension [36 that might be used to redefine the quantal 
system into a well behaved form that will render these interesting V function shape 
trajectories meaningful. In this way we actually compactify the a (or p) semiaxis, 
making these trajectories effectively periodic, in accordance with the periodicity of 
the V function. 

Any of the semiclassical trajectories p(t) with bounded amplitude from fig. (3) or 
their corresponding a(t) configurations in figs. (1) and (2) defines a vibrational mode 
of the O(N) vector. In such a vibrational mode, the unharmonic oscillator is described 
effectively by an harmonic oscillator that is coupled to a time dependent potential 
v(t) = + (3a{t) (Eq. Q2.4Q ). Each vibrational mode corresponds to a unique 
er(t) configuration found from Eqs. ( [2.16D or ( f2.32j ). On top of these vibrations, the 
O(N) vector rotates as well, giving rise to a host of rotational states that form a band 
superimposed on top of the vibrational mode. To leading order in 1/N, this rotational 
band collapses into a single energy vlaue as should be clear from Eq. ( p. 25 ). 

Each vibrational mode a(t) is an unstable saddle point of S e ff[a]. Indeed it can 
be shown straightforwardly that the second variation of S e ff[a] around each such 
a(t) configuration has an infinite number of negative directions, corresponding to the 
infinite number of decay modes of this excited mode as iV — > oo. Nevertheless, we 
must include these unstable saddle point configurations in calculating W P (T) in Eq. 
( |2.8| ) for finite T, since excited states of the oscillator do contribute to it. We will 
not pursue the case of the unharmonic oscillator any further in this paper, but rather 
turn to the two dimensional field theoretic case in the next section. 
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Ill The 0(N) Vector Model in Two Dimensions 



The action for the two dimensional O(N) model reads 



T oo 

S = J dt J dx 

-oo 



am 



&/3g 2 
AN 



(* 2 ) : 



(3.1) 



where $ is an N component scalar field transforming as an O(N) vector, m 2 is the 
bare "meson" mass and g 2 is the finite quartic coupling. a,/3 and 7 (that appears 
below) have the same meaning as in the previous section, hence (<9 M $) 2 = (<9o$) 2 — 
a(di<&) 2 . Rescaling 

$ 2 = Ncj) 2 



and introducing the auxiliary field a(x) as before, Eq. ([3.1D turns into 



S = N dt dx 

JO J -00 



1 , „ , , 2 am 



2 + a(3(a 2 - #a0 2 ) 



(3.2) 



(3.3) 



Note again, that by its equations of motion a = | (f) is a non-negative field (recall our 
convention g > 0). This positivity should persist in the quantal domain as well, since 
this equation of motion may be enforced equivalently as a constraint by introducing 
a Lagrange multiplier. 

In complete analogy with Eqs. (|2.6|)-([2T9|) we define the amplitude 



W P (T) = (0 O ; t = T|0 o ; t = 0) 
= J V pbc a e^/M 



(3.4) 



where </> has been integrated out yielding formally 

T 00 

[c] = a/3 J J dxa 2 Tr pf)C ln — □ — am 2 — 2af3ga 

-00 ^ 



(3.5) 



in which □ = d 2 — adj. 

As was mentioned in the introduction the ground state of the O(N) model is 
governed by the homogeneous constant o = o c configuration, where a c found from 
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ii i- • SS e ff[a 

by the extremum condition 



0, the so called "gap equation" of the 

model, and is a local minimum of S e ff[a}. This gap equation expressed in terms of 
the bare quantities m and g in Eq. ( |3.1| ) contains logarithmic divergences that can 
be swallowed into the bare m 2 parameter, replacing it by a renormalized finite m 2 R 
parameter, that depends on an arbitrary finite mass scale, while g 2 remains unchanged 



In general, this "gap equation" might have more than one solution a = o~ c , but only 
one such solution corresponds to the ground state of the theory. The latter is 0(N) 
invariant || as we are discussing a two dimensional field theory with a continuous 
symmetry []3"7[] . 

The local minimum a = a c of S e ff[a] supports small fluctuations of the fields 
in Eq. (|3.3|) around it giving rise to a low energy spectrum of "free mesons" with 
mass squared m R + 2(3ga c . There exist also bound states and resonances in the O(N) 
singlet (f)- 4> channel ||. This picture is intuitively clear in the case of positive quartic 
interactions (3 = 1, but it turns out to persist also in the (3 = — 1 case of negative 
quartic coupling |§, provided g 2 is not too large. What happens in the latter case 
is that despite the fact that for (3 = — 1 Eq. has no ground state, S e ff[a] still 

possesses a local minimum that support small oscillations of the fields in Eq. ( |3.1|) . To 
leading order in 1/N these fluctuations cannot spill over the adjacent local maximum 
of S e ff[a] towards infinite field values, provided g 2 does not exceed a certain maximal 
critical value. It was found in || that the a — a propagator is indeed free of tachyons 
even in the (3 = — 1 case, provided g 2 is smaller than that critical bound. 

We now turn to the central topic in this work, namely, the static a(x) extremal 
configurations of S e ff[a], and the bands of bound 4> states associated with them. Note 
that such static o~(x) configurations satisfy the periodic boundary condition along the 
time direction 



a(x, 0) = a(x, T) 
implied by Eq. ( |3.4|) automatically, for any value of T. 
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(3.6) 



Our approach to these static cr(x) configurations will be indirect. Namely, we 
assume that the extremal a(x,t) configuration we are looking for has a very slow 
time variation, such that [do, a] ~ will be a very good approximation in calculating 
the trace in Eq. ( |3.5| ). In this "adiabatic" approximation we effectively carry a 



smooth dimensional reduction of the extremum condition 5 f e j i ^ = from 1 + 1 to 
1 + space-time dimensions. Therefore, we use this "adiabatic" limit merely as an 
infra-red regulator of the t-integrations that occur in matrix elements in the trace of 
Eq. ( |3.5| ) which enables us to extract the static o~(x) = cro(x) piece. We do not know if 
the slowly time varying a(x, t) configurations we encounter in this approximation are 
relevant as genuine space-time dependent fields, or just artifacts of the approximation. 

Carrying the "adiabatic" approximation we expand the real field a in frequency 
modes 

oo 

a{x,t)= £ a n {x)e^ 1 (3.7) 

n=— oo 

where by virtue of Eq. (|3.6| ) and the fact that a* = a we have 



2ir 

uj n = — P{n) ■ P(-n) = -P(n) ; <(x) = a_ n {x) (3.8) 

in which P(n) is any monotonously increasing map from the integers to themselves. 
The precise form of P(n) will be immaterial to us, since we are interested only in the 
zero frequency mode ao(x) in Eq. (|3.7[). 

Using Eqs. and Q S effW\ in Eq. (pfe becomes 



SeffW] = a (3T J dx 



-oo 
oo 



Oq(x) + 2 o-_ n (x)cr r , 



x 



n=l 



1 r r i 

— dx ^2 (x,n\\n uj„ + ctd 2 x — am 2 — 2a(3ga n \x,n) (3.9) 



2 7 _ 

where {\x, n)} is a complete orthonormal basis of positions x and frequency modes 
u n and we have used the adiabatic approximation [8q, a] ~ 0. 

Varying Eq. (|3.9|) with respect to (k > 0) the extremum condition reads 
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for the k-th mode in Eq. (|3.7| ). Here we have used the relations (n\n) = 1 and a 2 = 1. 
Eq. ( 3.10| ) is of the general form of Eq. ( [2.11 ) and thus shows explicitly that our 
adiabatic approximation has reduced the 1 + 1 dimensional extremum condition into 
a + 1 dimensional problem, as was stated above. Indeed, Eq. ( |3.1U| ) for k = could 
have been obtained equivalently by ignoring completely any time dependence in the 
extremum condition 



5S eff [a] _ ^ 



2a(x, t) + - (x, t\ — - — — — — — — \x, t) 

7 —do + ~ am ~ 2apga 



0, (3.11) 



5a(x, t) 

namely discarding t, \t) and do in this equation from the very beginning. In order to 
maintain the correct scale dimensions of the remaining objects, one simply replaces 
g by g/T in Eq. ( |3.11|) , obtaining Eq. ( ft.lOQ again. 

The advantage in performing the adiabatic approximation (Eqs. ( ^3 . 7| ) - ( PTTOj ) ) 
rather than the simpler derivation of the extremum condition for <Jq{x) we have just 
describe is that it enables us to see explicitly what happens to non-zero frequency 
modes in Eqs. Q- 

At this point it is very important to note that we must exclude space independent 
solutions to Eq. ( |3.1U| ) (as far as k = is concerned). For such solutions Eq. ( 3.1U| ) 
is equivalent to an algebraic cubic equation similar to Eq. ( |2.12|) . However, the 
homogeneous a = a c vacuum condensate satisfies the "gap-equation" which is a 



transcendental equation H, |22j rather than an algebraic one. For this reason it is 



very clear that our resulting static (Tq{x) configuration will be very different from the 
one found in f2"2|| , which is only a mild distortion of the vacuum a c condensate. Thus, 
the band of bound states associated with our cr (x) complements those discussed in 
[Z2 . 

Since we are interested essentially in the zero frequency mode of Eq. ( |3.7| ), we are 
free to assume that a(x, t) is an even function of time, namely 
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cr_ fc (x) = a k {x) = a* k {x) (3.12) 

In this case, we see that Eq. ( 3.10 ) relates the diagonal resolvent R of the 
Schrodinger operator 



h = -d 2 + 2f3ga k (x) 
at energy parameter oluj\ — m 2 , to its potential term a k (x): 



, 2cryT 

R(x) = a k (x) 
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(3.13) 



(3.14) 



in complete analogy with Eq. ( [2.13 



Following the same steps as in the previous section, we find from the Gelfand-Dikii 
equation (Eq. A. 12) that a k must satisfy the equation 



2a k a' k ' - {a' k f - 4 (2(3ga k + m 2 - au 2 k ) a 



2\ Jl 
k 



a 



(27) 



(3.15) 



where we have used the relations a 2 = 1 and 7 2 = —a. Note from Eq. ( |3 . 1 U| ) that Eq. 
( |3.15|) holds for the zero frequency mode a (x) independently of whether the adiabatic 
extension o~(x,t) of the static configuration has any definite parity properties under 
t — > — t or not. 



Naturally, Eq. (|3.16|) coincides with Eq. (|2.14|) in the case of euclidean signature 
a = — 1 upon the replacements m 2 — - > m 2 + uj 2 , g — > g/T. 

As in the previous section, first integration of Eq. (|3.15|) is done by substituting 



Ou (X 



fk(cr k ) 



(3.16) 



leading to 



f, 



' da k 
dx 



Pgo-l + ( 



m — auo h a 



o_k_ 



(3.17) 
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where a is an integration constant. Note that unlike the + 1 dimensional model 

discussed in the previous section, in the field theoretic case we lack a simple physical 
interpretation of Eq. ( |3.17| ) as an equation of motion in one dimension allowing us to 
discard "radial" trajectories (in euclidean signature) that hit the origin a = 0. 



Simple analysis of Eq. (|3.17|) for very high frequencies C0k, reveals that the corre- 



sponding modes ak(x) with bounded amplitudes are highly suppressed. Namely, one 

a 2 



2 

finds that is of the order —. 9 )2 -C 1. Therefore, if we choose P(n) in Eq. 



such that already 2tc -P(l) ^> 1 (e.g. P(n) = ne'" +1 < l ), then for large enough T 
and generic values of a all non zero frequency modes will be highly suppressed in the 
resulting extremum condition, assuring self consistency of the adiabatic approxima- 
tion. As cro(x) is decoupled from non zero frequency modes within the framework of 
the approximation made, we are free to make such choices of the P(n)'s. 
Concentrating on the zero-frequency mode, Eq. ( |3.17| ) becomes 



da 

dx 



2 



= f{a) = 4 f3ga A + mV + a ( ) ( ^ - 1 ) (3.18 



V2T 



a 



a 



where we have dropped the sub-index k = everywhere. Note the explicit T de- 
pendence of the right hand side of Eq. ( |3.18|) . This implies that o"o(x), which is our 
time-independent saddle point configuration, will have a peculiar dependence on T. 

However, this paradox will be resolved by the fact that solutions to Eq. ( |3.18| ) 
in which the T dependence is "important" (in a sense that will become clear when 
we consider these solutions) always give rise to infinite S e ff[a] values, and are thus 
highly suppressed in the path integral ( Eq. ( |3.4| ), unless we put the field theory in a 
finite "spatial" box of length L. In the latter case, these T dependent cr (x) solutions 
are not suppressed. In such cases having explicit T dependence is consistent with 
having explicit L dependence. We will not discuss the interesting finite size behavior 
of this field theory and its corresponding thermodynamic behaviour in this paper. 

Thus, if <7i, o"2 and are the three roots of the cubic equation / 2 (0) = 0, we have 
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/ 2 (0) = -«(^) 2 (3.19a) 
/ 2 (±oo) = ±(3oo (3.19b) 
oio^ = — (3.19c) 

(71 + ^2 + ^3 = • (3.19d) 

9 

Here we also adopt the convention that if o"i, a 2 and 03 are all real, then 

o"i < 02 < o"3. Unlike Eq. ( |2.16| ), a is not an overall factor in ( |3.18| ) and we shall 



have to analyze all four possible (a, (3) combinations separately. We discuss below 
only cases in which a% , 02 and 03 are all real. 



Clearly, then, ( |3.19| c) implies that each a, (3 combination gives rise to two possi- 
bilities, according to the arithmetic signs of o"i,cr 2 and <7 3 , yielding eight possibilities 
in all. These eight (Jq(x) configurations are depicted schematically in figs. (4a)- (4d). 
However, two of these (figs. (4a.i) and (4b. i) are explicitly ruled out by (|3.19| d), since 
it implies that when o"i,o"2 and cr 3 have all the same sign, the latter must be —f3. 
In the remaining configurations of fig. (4), the physically allowed domains for a are 
segments on the semiaxis a > along which f 2 (a) > 0. 

As in the previous section, it is clear from ( |3.18| ) and figs. (4a-d) that all static a 
configurations are of the form u{x) = ip(x — Xq) where xq is an integration parameter 
and ip is either (essentially) a Weierstrass V function (unbounded a configurations) 
or a simple rational function in terms of Jacobi functions, e.g. the sn or cn functions 
(bounded a configurations). All these a configurations are periodic in x. In euclidean 
space-time signature we have again the a configurations that hit the endpoint o = 0. 
Unlike the + 1 dimensional case discussed in the previous section, we have no a 
priori reason to discard such solutions, but note that as functions of x, their first 
derivative suffers a finite jump discontinuity whenever a vanishes. This fact makes 
such cr(x) configurations rather suspicious. 
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Fig. 4c a = - (3 = - 1 




Fig.4d a = -p = -l 
FIGURE 4 
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We will discuss explicit solutions of Eq. ( |3.18| ) in the last part of this section. At 
this point we make a small digression to recall the way the spectrum of "mesonic" 
bound states is calculated, given an extremal static a(x) configuration (i.e., a specific 
solution to Eq. (|3.18|) in our case). In this digression we follow [20 , p2| very closely. 
Given a static a configuration and ignoring the gaussian fluctuations of a around it, 
dynamics of the "meson" (f) field in this background is governed by the action in Eq. 
(|3.3|) upon substituting into it the particular a(x) configuration we are discussing. 
Clearly then the fluctuations of (f) satisfy 

[□ + am 2 + 2af3ga(x)}(j) a = 0; a = l,...N (3.20) 

with the periodic boundary conditions 

(f> a (t + T,x) = (f> a (t,x). (3.21) 

Upon a simple separation of variables (fi a (x,t) = e lXt (j> a {x) Eq. (|3.20j ) turns into the 
eigenvalue problem of the time independent Schrodinger equation^ 

-d 2 x +m 2 + 2(3ga(x)] <j> a (x) = A 2 [or]0 a (x) (3.22) 

where the functional dependence of the eigenvalue A 2 on the Schrodinger potential 
cr(x) has been written explicitly. 

Following Eq. ( |3.22[) we expand </> a (x) as 

a (x,t)=^^(t)^(x) (3.23) 

i 

where {ipi{x)} are the complete orthonormal set of eigenfunctions of the Schrodinger 
hamiltonian in Eq. ( |3.22| ) with eigenvalues {A 2 }, and A % a {t) are arbitrary amplitudes 
to be integrated over in the functional integral. The latter are subjected to the 
periodic boundary condition A l a (t + T) = A l a {t). 

Using Eqs. ( |3.20| )-( |3~2~3D the (f) dependent part of Eq. ( p^3|) becomes 

N T - N 



5 = tMEE [^Kf - *\ 2 {A\Y] (3.24) 

q a=l i 



22 We have omitted an overall factor a from one of the sides of Eq. (3.21). This will show up 
later. 



31 



which is nothing but the action iV identical (infinite) sets of harmonic oscillators 
oscillating at frequencies {A,}. 

The natural place to look for bound state energies is the partition function as- 
sociated with Eq. ( |3.24|) , i.e. with the particular static a(x) configuration we are 
considering. 

Using well known results |20], ^ to integrate over the A l a {t) in Eq. ( |3.24|) , the 
desired partition function associated with a(x) 



Z[a] = Tre 11 ^™ = J V pbc (j> W P (T; [a]) 



(3.25) 



(where Wp(T; [a]) is the contribution of our particular a(x) configuration to Eq. 
(|3.4|) ) may be expressed straightforwardly as 



Z\a] 



n i _ e \ lT h 



N a/3~/NT j a 2 dx 



■ e 



NT 



EC({n})e 

{n} 



J2j n 3 X 



T 

7 ^3 



(3.26) 



The sum in Eq. ( 3.26 ) runs over all ordered sets {n} of non negative integers 



and 



c(w)= n 



(N + m-i)\ 



By definition we have (Eq. ( |3.25| )) 



Z[a] = £C({n})e E «"» T ^ 

{n} 



whereQ 



E({n}) 



/oo AT 
a 2 dx + —J2 X 
-oo 2 



+ E n 3 X J 

n,je{n} 



(3.27) 



(3.28) 



(3.29) 



is the energy of the state corresponding to the set {n} ^ and C({n}) is its degen- 
eracy. 



23 We use the relations 7 2 = —a, a 2 = 1. 

24 That is, inserting "mesons" into energy level Xk , k = 1, 2, . 



25 Clearly E({n}) corresponds to a " meso nic" bound state only if all quantum numbers that 
correspond to scattering states in Eq. ( 3.22 ) vanish. 
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Therefore each cr(x) saddle point configuration gives rise to a band of bound 
states (as well as scattering states) of <fi "mesons" with a very rich O(N) structure. 
Generically, the degeneracy factors of energy states in the band are very large and 
each E({n}) multiplet forms a highly reducible representation of O(N). The lowest 
state in a band has no "mesons" in it at all (rij = for all z's) and its energy is given 
by 

7 N 

E [a] = -f3N J a 2 dx + y E X M ( 3 - 30 ) 

— oo % 

where the first term is the classical energy of the o[x) field and the other term is the 
zero point contribution of the oscillators in Eq. fl3-24|) . 

In general both terms in Eq. ( |3.30| ) are divergent and must be regularized. In 
certain cases where the period of a[x) as a function of x becomes infinite regular- 
ization of Eq. ( |3.30| ) may be achieved simply by subtracting from it the zero point 



contribution of oscillators in the background of the vacuum configuration a = o~ c , 
namely the quantity J2i \[cc) |^0, E^j-H We elaborate on such cases below. In other 
cases E Q [a] might be inherently divergent, suppressing contributions from such o~(x) 
configurations enormously relative to the vacuum. Such cases include the Weierstrass 
V function configurations as well as finite amplitude <r(x)'s with a finite period along 
the x axis. However, the latter may become important if we put the field theory into 
a spatial "box" of finite length L, that is an integral multiple of period of a(x). As 
L — > oo these modes are suppressed as was stated above. 

In [2D, |22| Eqs. ( |3.29| ) and ( |3.30| ) were the starting point in a calculation determin- 



ing extremal a(x) configurations. In these papers one extremises the regularized form 
of Eq. ( |3.29| ) with respect to a(x). In order to account for the functional dependence 
of the Aj on a(x) one has to invoke inverse scattering techniques. In these papers, 
the <j(x) configurations that resulted from the inverse scattering analysis were generi- 
cally very mild distortions of the vacuum condensation a = a c obtained from the gap 
equation, expressed in terms of reflectionless Schrodinger potentials. Since a = a c 
is a local minimum of the appropriate S e ff[a], it is clear from the self consistency of 
26 Note in this respect that our 2ga c corresponds to [x~ X)/N in [0. 
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these analyses that such cr(x) configurations must be also local minima of S e /f[a}, 
since they are just the exact result (as iV — > oo) for the back-reaction of either bosons 
or fermions on the vacuum, which is in turn the local minimum of S e ff[a] trapping 
these particles. 

As was mentioned above briefly, our static extremal cr(x) configurations are very 
different from those found in p0|, [22J. Moreover, contrary to the latter, they are 
generically not local minima of S e ff[a] but rather saddle points of this object such 
that the second variation of S e ff[o~] around them is not positive definite. We have 
seen in the previous section that despite this fact the contribution of these saddle 
points to the path integral representation of the partition function Tr e HT ^ is impor- 
tant, because they correspond to (collective) excited states of the quantized theory. 
Clearly, only the extremal configuration that describes the ground state of the quan- 
tum theory must be a local minimum of S e ff[a]. We may strengthen our conclusion 
of the importance of these o~(x) configurations by making an analogy with a simpler 
elementary example. Namely, consider a variational calculation of the spectrum of a 
Schrodinger hamiltonian (for simplicity we consider the one dimensional case) 

h = -d 2 x + V(x) . (3.31) 

The energy functional is 

b 

'ifj\ iJr,E]= J ip\h - E)ip dx + E (3.32) 



S 



where E is a Lagrange multiplier enforcing the normalization condition (^l^) = 1 
and Neumann boundary conditions are assumed at the end points. 
The extremum condition s ^f^ = reads 

(h — E)ip = ; ^'(a) = ^'( b ) = • (3.33) 

Obviously, Eq. ( |3.33| ) is nothing but the Schrodinger equation of h, and thus has 
an infinite number of solutions {ip n (x),E n } (n = 0,1,2...) - the spectrum of h. 
Clearly, only one extremum is an (absolute) minimum of 5* - this is the ground state 
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{i[!q(x), E }. The other extrema corresponding to excited states are only saddle 
points of H. Indeed, the n-th state, ip n (x), has exactly n unstable directions in its 
saddle point (recall we are considering the one dimensional case where no degeneracy 
occurs), corresponding to the n states below it. This is clear from the second variation 
of S around {i[) n (x), E n } 

oo 

S [V>£ + + <ty, E n ] = J2 \C m \ 2 (E m - E n ) + E n (3.34) 

m=0 

where Sip = J2m=o C m ijj m {x) and we have used the orthonormality of the {ip n }- These 
saddle points are important as they correspond to all excited states of the spectrum 
and must be taken into account. 

Having established the importance of our static saddle point configurations we now 
turn to the explicit calculation of these. Out of eight cases depicted in fig. (4) only 
two yield a solution to Eq. ( |3.18| ) that can produce a finite S e ff[a] value. Therefore, 
these are the only cases relevant for our discussion of bound states and resonances. 
These cases are described by figs. (4b. ii) and (4d.i), which have the same behaviour 
around a 2 and 03. As will be shown below, the finite S e ff[a] value for these cases is 

obtained in the limit a 2 — > o~\-\ > 0. Two of the other cases in fig. (4) have already 

been ruled out by Eq. (|3.19| d) and the remaining four give rise to infinite 5 e //[cr]'s or 
contain o~(x) configurations in euclidean time that hit the a = boundary point. 

In both cases of interest to us here the quartic interaction is unbounded from 
below {(3 = —1), but as was discussed at the beginning of this section, the field 
theory is perfectly defined to leading order in 1/N in terms of the fluctuations around 
the vacuum of the theory which is a local minimum of S e ff[o~]. Though our resulting 
cr(x) configurations are quite distinct from the vacuum condensation itself, they are 
of a finite amplitude, and are therefore protected against being driven to infinity by 
the bottomless quartic interaction. The solution of Eq. ( p.l8|) for both these cases 
read 



a(x) = a 2 + (cr 3 - a 2 )cn 2 



\fgio3 - o"i) (x - x ) \n 



(3.35) 
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where cn is a Jacobi elliptic function with parameter [|4], |33 

/x = • (3.36) 

For generic values of a\, 02 and 03 we have < \i < 1 and a(x) has a finite real 
period L along the x axis given by 



2K(n) 



(3.37) 



where K(n) is a complete elliptic integral of the first kind. In this case Eq. ( |3.35|) 
yields obviously an infinite S e ff[a] value, unless the field theory is defined in a "box" 
of finite volume. However, for /i — > 1 — , L in Eq. ( 3.37|) diverges as | In and cr(x) 
degenerates into 

<r(x) w o- 2 + 32 sech 2 [y^E (x - x )] +0(1-/1) (3.38) 
n-+i- 

where = 03 — 0^; k = 1, 2. Eq. (|3.38| ) leads evidently to a finite S e ff[a) value. 

The desired limit \i — > 1— is obtained as <ti — > o"2— . In the Minkowsky space 
theory (fig. (4b. ii)) this limit is possible only as both <j\ and o"2 vanish, since they 
have opposite signs. This is possible only when {^§f\ in Eq. ( |3.18| ) becomes very 
small. In such a case we have from Eq. ( |3.18| ) 



m 

a 3 = — + ( 

9 

0"2 = —G\ = 

Using Eq. ([U|), Eq. (|3T38| ) becomes (for ft = 1) 

gcr(x) = m 2 sech 2 [m(x — xq)} . 





m 


TV 


9 


(f) 


w 





(3.39) 



(3.40) 



Strictly speaking the euclidean space theory approaches this limit through complex 
0- 2 = o"i values, and fig. (4d.i) does not hold in this case, but the end result, Eqs. 
( |3.39| ) and ( |3.40p are the same, as is obvious from Eq. ( |3.18j ). 

Substituting Eq. ( ^401) into Eq. (|3T2p (for (3 = -1, a = ±1) we find that (f> is 
coupled to the one dimensional potential V(x) given by 



V[x) = m 2 — 2ga(x) 



m 



2m 2 sech 2 [m(x — xq)] . 



(3.41) 
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V(x) is a reflectionless Schrodinger potential with a single bound state ipo(x) at zero 
energy A = f40fl , where 



^o( x ) = v/ secn i m ( x ~ x o)\ ■ {3A2) 

In addition to if>o(x), V(x) has a continuum of scattering states 

ip q (x) = {iq - m tanh [m(x - x )]} e iqx (3.43) 

with eigenvalues = m? + q 2 in which the reflectionless nature of V(x) is explicit. 

In this respect, our u{x) configuration is reminiscent of the Callan-Coleman-Gross- 
Zee kink in the Gross- Neveu model ||2U|| . As in that kink solution, the profile of a(x) is 
independent of the number of particles (<fi "mesons") trapped in it. This is contrary to 



the behaviour of cr(x) configurations found in p0[ by inverse scattering techniques, 
that are all small distortions of the vacuum configuration a = a c . We would like to 
stress at this point that unlike the Callan-Coleman-Gross-Zee kink, our configuration 



(Eq. ( |3.41|) ) does not connect degenerate vacua, since the O(N) model lacks such 



structures. 

As far as "mesonic" bound states are concerned, we have only one relevant quantum 
number uq in Eqs. ( |3.26|) - (|3.29| ), corresponding to Eq. ( |3.42j ). In this case the 



degeneracy factor in Eq. (|3.27|) reads 

CM = l^LZill ,3.44) 
(N - 1)! n \ 

which implies that the dimension of the multiplet with eigenvalue E(uq) is that of a 
symmetric tensor of rank n which is a highly reducible representation of O(N), the 
irreducible components of which are all traceless symmetric tensors of lower ranks n' 
such that n — n' Q = (mod 2). 

Up to this point we have not specified the value of m in any of our equations. 
Explicit comparison with Eq. (4.38) of |22j as well as the analogy with the Callan- 
Coleman-Gross-Zee kink indicate that we must identify m as the finite (renormalized) 
mass of the (f> "mesons" obtained from the "gap-equation" 0. Moreover, from Eq. 
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(4.14) of |22[] it seems that our "kink-like" cr(x) configuration corresponds to the 
parameter 9 in that equation having the value 9 = |. Ref. p2] considered only 
the range < 9 < |, discarding 9 — |. The reason for this was that only in that 
region did E(n Q , 9) obtained extrema which were actually local minima of E(n Q , 9) as 
a function of 9. But since we are discussing excited states of the field theory, it is clear 
that these minima prevail only in a subspace of the total Hilbert space of states that 
is orthogonal to all lower energy states (recall Eqs. (|3.31| )- (|3.34j )). Therefore there 
is nothing wrong with our o~(x) configurations lying outside the domain considered 
by It has been shown in [^2| that in the case of stable quartic interactions 

(j3 = +1) no saddle points a(x) configurations that are minor distortions of the 
vacuum configuration a = a c could be found. We complement this conclusion by 
proving that in this case "big" kink like c{x) configurations do not occur as well. 

We close this section by actually calculating the spectrum given by Eq. ( |3.29|) 
and (|3.30|) . The regularized lowest energy state in the bound is given by 

E^[a] = E [a] - E [m] , (3.45) 

Namely, following Eq. (|3.3(J| ) (for (3 = — 1) 
1 °° 1 

- E***[a] = J a 2 dx+-J2 (A,W - A,[m 2 ]) . (3.46) 

Note that V(x) in Eq. ( |3.41| ) is isospectral to the "potential" V(x) = m 2 , except for 
the single zero energy bound state of V(x) ||40|| . Therefore, all terms in the sum over 
Aj in Eq. ( |3.46| ) cancel identically except for Ao, which is zero anyway. Thus, only 
the integral in Eq. (|3.46|) contributes to E^ eg , yielding 

*fH = ^. (3.47) 

Since the only discrete eigenvalue of V(x) in Eq. ( |3.41| ) is zero, Eq. (|3.47| ) yields 

the common value of all bound state energies 

, s 4Nm 3 
E(n ) = (3.48) 

independently of no- In this case the whole band degenerates into a single energy 
level. 
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Conclusion 



In this work we have shown how simple properties of a one dimensional Schrodinger 
hamiltonian can be used to obtain time or space dependent saddle point configurations 
contributing to the path integral in the large N limit. As a specific model we have 
analyzed the 0{N) vector model in + 1 and in 1 + 1 space-time dimensions. 

In the quantum mechanical case our results are equivalent completely to those 
obtained directly from the radial Schrodinger equation in the large N limit. The 
time dependent o~(t) configurations we have found are nothing but the semiclassical 
radial trajectories of that equation. In this respect our calculation is to the radial 
Schrodinger equation exactly the same as what the instanton calculus is to JWKB 
calculations of tunnelling effects. 

In the two dimensional field theoretic case our method has been applied to find 
static saddle point configurations. In this case, the static case was approached through 
an "adiabatic" approximation that discarded all time dependence. This led to static 
saddle points that were quite distinct from the homogeneous vacuum one and from 
other static configurations obtained as minor distortions of the vacuum by use of 
inverse scattering techniques. Many of our novel static saddle points are important 
only in finite volume cases. We did not elaborate on these in this paper. Only 
one homogeneous saddle point configuration turned out to be relevant in the case of 
infinite volume. The latter gave rise to "mesonic" bound states at zero binding energy 
and is reminiscent of the Callan-Coleman-Gross-Zee kink found in the Gross-Neveu 
model. 

All extremal o~(x) configuration we have found in the O(N) vector model turn out 
to be saddle points of S e ff[a] rather than local minima of that object. This, however 
should not be of any surprise, since our saddle points correspond to excited states of 
the quantal system. 

Our method could be applied straightforwardly to other two dimensional field 
theoretic models. 
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Appendix 



The differential equation obeyed by the diagonal resolvent of a Schrodinger 
operator in one dimension 

Consider the spectral problem for the one dimensional Schrodinger operator 

h = -d 2 x + U(x) (A.l) 

on the segment a < x < b (where any of the end points may be at infinity) with 
boundary conditions 

aip(z) + Pif>'(z) = 0; z = a, b (A.2) 

on the wave functions. 

The eigenfunctions ip n {x) and the eigenvalues E n of h are therefore the solutions 

of 

-dl + U(x)] ip(x) = E ijj(x) (A.3) 

subjected to the boundary condition of Eq. (A.2). Let E be any real number, and let 
ip a {x) (ipb(x)) be the solution to Eq. (A.3) that satisfy the boundary condition Eq. 
(A.2) only at x — a (x — b), but not at the other boundary point. 
The Wronskian of ipa and ipb 

W ab = V«(aOV£(*) - 1>a(x)Mz) = C(E) (A.4) 

is x independent, and clearly vanishes (as a function of E) only on the spectrum 
of h For a generic E, the Green function corresponding to Eq. (A.3) is the 
symmetric function. 

G E (x,y) = [e(x - y)i[) a (y)i) b (x) + Q(y - x)i) a (x)i) b (y)} (A.5) 



27 Indeed, if a and b are finite, the spectrum of h is non-degenerate. If either a and/or b become 
infinite, the continuum part of the spectrum (if exists) is only two-fold degenerate. Explicitly, if E 
is an eigenvalue, ip a must satisfy the boundary condition at x = b as well, leading to C{E) = 0. 
Conversely, if C(E) = 0, ip a {x) satisfies the boundary condition at x = b, hence E is in the spectrum. 
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where Q(x) is the step function. 

It is easy to check that G E (x,y) satisfies its defining equation 



-d 2 x + U(x) G E (x, y) = -dl + U(y) G E {x, y) = S(x - y), (A.6) 



and therefore 



G E (x,y) = (x\ 



1 



-d 2 + U-E 



The diagonal resolvent of h at energy E 



Re(x) = (x\ 



1 



—d 2 + U — E 



\x) 



is obtained from as 



R E {x) = hm - [G E (x, x + e) + G E (x + e, x) = — — — — 

e-*0+ z G [E) 



Using Eq. (A. 3) we find 



R'Ux) = 2 



(U-E) + 



R 



leading to 



R'l - 2R E R" E = -AR 2 E (U — E) + I -7—7- R 2 . 



Using eqs. (A. 4) and (A. 9) we finally arrive at 



(A.7) 



(A.8) 



(A.9) 



(A.10) 



(A.ll) 



2R E R" E + R' 2 + 4R 2 E (U -E) = l 



(A.12) 



which is nothing but the "Gelfand-Dikii" equation [23|, alluded to in sections II and 
III. 
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